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Abstract 

Context. The detection of reflected light from an exoplanet is a difficult technical challenge at optical wavelengths. Even though this 
signal is expected to replicate the stellar signal, not only is it several orders of magnitude fainter, but it is also hidden among the stellar 
noise. 

Aims. We apply a variant of the cross-correlation technique to HARPS observations of 51 Peg to detect the reflected signal from 
planet 51 Peg h. 

Methods. Our method makes use of the cross-correlation function (CCF) of a binary mask with high-resolution spectra to amplify 
the minute planetary signal that is present in the spectra by a factor proportional to the number of spectral lines when performing 
the cross correlation. The resulting cross-correlation functions are then normalized by a stellar template to remove the stellar signal. 
Carefully selected sections of the resulting normalized CCFs are stacked to increase the planetary signal further. The recovered signal 
allows probing several of the planetary properties, including its real mass and albedo. 

Results. We detect evidence for the reflected signal from planet 51 Peg b at a significance of 3 -(T'„oi 5 j. The detection of the signal 
permits us to infer a real mass of 0.46);|5 qj Mj„p (assuming a stellar mass of 1.04 Mq) for the planet and an orbital inclination of 80);]® 
degrees. The analysis of the data also allows us to infer a tentative value for the (radius-dependent) geometric albedo of the planet. The 
results suggest that 51Peg h may be an inflated hot Jupiter with a high albedo (e.g., an albedo of 0.5 yields a radius of 1.9 ± 0.3 Rj^p 
for a signal amplitude of 6.0 ± 0.4 x 10“^). 

Conclusions. We confirm that the method we perfected can be used to retrieve an exoplanet’s reflected signal, even with current 
observing facilities. The advent of next generation of instruments (e.g. ESPRESSO@VLT-ESO) and observing facilities (e.g. a new 
generation of ELT telescopes) will yield new opportunities for this type of technique to probe deeper into exoplanets and their 
atmospheres. 
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1. Introduction 

Since the discovery in 1995 of a planet orbiting 51 Peg (Mayor 
& Queloz 1995), over 1800 planets in around 1100 planetary 
systems have been found (Schneider et al. 2011)*: this number 
increases steadily. Furthermore, close to 470 multiple planetary 
systems have been detected, some of which are highly complex 
(e.g., Borucki et al. 2013). 


* Based on observations made with ESO Telescopes at the La 
Silla Paranal Observatory under programme ID 09LC-0271 (with the 
HARPS spectrograph at the ESO 3.6-m telescope). 

** The radial velocity data for the HARPS observations of 51 Pegasi 
are only available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
bin/qcat?J/A-tA/ 

* http://exoplanet.eu/ 


The search for exoplanets has been following two different, 
but complementary, paths: the detection of exoplanets with in¬ 
creasingly lower masses, and the characterization of these exo¬ 
planets and their atmospheres. On the detection side, the radial 
velocity (e.g., Lovis & Fischer 2010; Bonfils et al. 2013) and 
transit methods (e.g., Borucki et al. 2010; Winn 2010) have been 
the most prolific. One of the most important resul?ts of planet de¬ 
tection surveys is the ubiquity of planets around solar-type stars 
(e.g., Howard et al. 2010). 

On the characterization side, the current frontier of exoplanet 
characterization has been pushed toward the study of exoplanet 
atmospheres, both from a composition and from a dynamics 
point of view. To overcome this difficult challenge, several indi¬ 
rect techniques have been developed. Transmission spectroscopy 
relies on observing the host star spectrum as it is filtered by 
a planet atmosphere during a transit (e.g., Charbonneau et al. 
2002; Knutson et al. 2014). Occultation photometry and spec¬ 
troscopy measure the occultation (secondary eclipse) depth of 
the star-Hplanet light curve at different wavelengths to derive the 
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planet’s thermal (e.g., Deming et al. 2005; Snellen et al. 2010; 
Stevenson et al. 2014) and reflected (Alonso et al. 2009; Rodler 
et al. 2013) signatures. The detection of the exoplanet emission 
spectra was also possible through high-resolution spectroscopy 
(Brogi et al. 2012; Rodler et al. 2012; de Kok et al. 2014). The 
measurement of phase variations relies on the detection of the 
flux variation along the planet’s orbit as it alternately presents its 
day and night hemisphere to us (e.g., Knutson et al. 2009; Kane 
et al. 2011). These techniques represent the current front line of 
exoplanet characterization and are limited only by flux measure¬ 
ment precision that they impose as a result of the low planet-star 
flux ratio (e.g., in the most favorable cases, for a Jupiter sized 
planet with a 3 day period orbit T'pianet/J^ 5 »ar ~ 10^“* in the visi¬ 
ble and Fp\^netlFStar ~ 10“^ in the IR). 

Charbonneau et al. (1999) attempted to recover the reflected 
spectrum of the giant planet orbiting t Boo, which paved the way 
for detecting reflected light in the optical. To do so, they per¬ 
formed a evaluation of simulated star-i-planet spectra against 
high-resolution observations obtained with the HIRES spec¬ 
trometer at the Keck observatory. Although unable to detect the 
reflected signal from planet t Boo b, they were able to set an 
upper limit on the planet’s maximum planet-to-star flux ratio of 
about 5 X 10“^. In the same year. Collier Cameron et al. (1999) 
made an attempt at detecting the reflected signal of the same 
planet using a least-squares deconvolution technique, also with¬ 
out conclusive results. More recently, other attempts with alter¬ 
native methods have been made (e.g., Leigh et al. 2003; Rodler 
et al. 2010), albeit all results have been inconclusive about a 
definitive detection of a reflected exoplanet signal. Nonetheless, 
all these attempts are of great scientific values because they al¬ 
low establishing upper limits on the planet-to-star flux ratio. 

More recently, researchers were able to collect measure¬ 
ments of the albedo of several exoplanets in an effort to constrain 
current atmosphere models (e.g., Cowan (& Agol 2011; Demory 
2014) and infer their composition (e.g., the presence of clouds 
in the atmosphere of HD 189733 b as done by Barstow et al. 
2014). An interesting result is the observation of HD 189733 
b, where researchers were able to infer the planet’s color from 
albedo measurements (Evans et al. 2013). Nonetheless, several 
of these results are still subject to some discussion because it 
has been shown that the blue excess in the planet’s transmission 
spectrum might also be the result of stellar activity (Oshagh et al. 
2014). 

The detection of sodium absorption in the transmission spec¬ 
trum of HD209458 b (Charbonneau et al. 2002) paved the way 
for the detection of spectral features in exoplanet atmospheres. 
As new facilities were developed and analysis techniques per¬ 
fected with time, other elements or even molecules were de¬ 
tected (e.g., Redfield et al. 2008; Snellen et al. 2010; Birkby et al. 
2013; Madhusudhan et al. 2014; Brogi et al. 2014; Snellen et al. 
2014). 

In this paper we apply the technique described in Martins 
et al. (2013) to HARPS observations of 51 Peg to try to re¬ 
trieve the reflected spectrum of its planetary companion. 51 Peg 
(HD217014) is a solar-type star, slightly more massive than our 
Sun {Msi Peg!M q ss 1.04, Santos et al. 2013), at a distance of 
approximately 15.6 parsec from us (van Leeuwen 2007). With a 
minimum mass slightly under half the mass of Jupiter and an or¬ 
bital period slightly longer than days, 51 Peg b is the prototype of 
a hot Jupiter, giant gas planets in close orbits (Mayor & Queloz 
1995). The combination of the host brightness (VMag - 5.46, 
van Belle & von Braun 2009), the giant planet’s large dimen¬ 
sions, and the short-period planetary orbit yield a relatively high 
planet-star flux ratio and make this planetary system an excel¬ 


lent candidate for testing the method presented in Martins et al. 
(2013). 

In Sect. 2 we describe the principle behind our method. In 
Sect. 3 we describe the method and its application to our data. 
The results are presented in Sect. 4 and are discussed in Sect. 5. 
We conclude in Sect. 6. 


2. Principle behind the method 

Martins et al. (2013) showed that the cross-correlation function 
(hereafter CCE) can be used to mathematically enhance the S/N 
of our observations to a level where the extremely low S/N plan¬ 
etary signal can be recovered. The CCE of a spectrum with a 
binary mask (Baranne et al. 1996) has been extensively tested in 
detecting exoplanets with the radial velocities method. Briefly, 
this technique corresponds to mapping the degree of similarity 
between the stellar spectrum and a binary mask (representing 
the stellar type), which increases the S/N of the data by a factor 
proportional to the square root of the number of spectral lines 
identified in the mask. 

As discussed in Martins et al. (2013), we expect the planet’s 
signature to replicate the stellar signal, scaled down as a result 
of geometric (e.g., planet size) and atmospheric (e.g., albedo) 
factors. The planetary albedo measures the fraction of incident 
light that is reflected by the planet atmosphere. Several albedo 
definitions exist (e.g., Marley et al. 1999; Collier Cameron et al. 
2002; Seager 2010), but in our study we only considered the 
geometric albedo Ag, which is defined as the ratio of a planet’s 
flux measured at opposition (a - 0) by the flux of a Lambertian 
disk at the same distance and with the same radius as the planet. 
This allows us to define the planet/star flux ratio as 


Fplanet(.^^ 
F Star 


'■ g{a) 


R 


planet 


( 1 ) 


where Ag is the planet’s geometric albedo, a the orbital phase, 
g{a) the phase function, and Rpianet and a the planetary radius 
and orbital semi major axis.. 

Eor 51 Peg b, assuming an albedo of 0.3 and a planetary 
radius of 1.2 Rjup, Eq. 1 will yield a maximum planet-star flux 
ratio of ss 2 X 10“^. To detect the planet’s reflected signal under 
these conditions, we consequently need data with a combined 
S/N of at least 165 000 for a 3 cr„oise detection. A typical G2 
stellar CCE-mask used in HARPS has over 4000 spectral lines. 
This means that in principle a spectrum with a S/N of about 2600 
will contain enough information to allow us to build a CCE with 
the necessary S/N to detect the light spectrum of 51 Peg reflected 
on its planet (at 3-cr level). A lower S/N will suffice if the albedo 
or the planetary radius are higher (according to Eq. 1). 


3. Method 

3.1. Data 

Our data were collected with the HARPS spectrograph at ESO’s 
3.6 m telescope at La Silla-Paranal Observatory as part of the 
ESO program 091.C-0271. It consists of 90 spectra observed 
in seven different nights (2013-06-08, 2013-06-25, 2013-08-02, 
2013-08-04, 2013-09-05, 2013-09-09, and 2013-09-30), which 
amounts to about 12.5 h of observing time. These observations 
were split into several carefully selected time windows in which 
the planet could be observed close to superior conjunction (i.e., 
when the day side of the planet faces us) to maximize the planet’s 
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Table 1. Available data for the individual nights. 


Night 

Number of 
Spectra 

Total exposure 
[seconds] 

S/N range 

2013-06-08 

3 

1360 

243 - 296 

2013-06-25 

10 

5260 

273 - 351 

2013-08-02 

2 

1092 

145 - 151 

2013-08-04 

20 

10000 

215 - 311 

2013-09-05 

4 

2400 

191 - 248 

2013-09-09 

13 

7810 

122 - 265 

2013-09-30 

39 

17100 

179 - 388 


flux (maximum phase). These time windows were computed 
from the ephemeris provided by Butler et al. (2006). 

The obtained spectra have a S/N on the 50* order (~5560A) 
that varies between 122 and 388. The spectra cover the wave¬ 
length range from about 3781A to 6910A. More detailed infor¬ 
mation can be found in Table 1. 

Despite the brightness of the target, some cloudy nights de¬ 
creased the expected S/N (e.g., S/N~150 after an exposure of 
600s on 2013-08-02, while on 2013-09-30 we managed to ob¬ 
tain ~390 after 450s). 

3.2. Data reduction 

To reduce the data and create the CCFs for each spectrum, the 
HARPS DRS (Mayor et al. 2003) was used. The data were re¬ 
duced using the default settings and were then fed to the CCF 
calculation recipe, used with a weighted G2 mask (Pepe et al. 
2002 ). 

Selecting an optimized CCF computation step was of critical 
importance. During the detection process, the CCFs need to be 
shifted in radial velocity (see below), which implies an interpola¬ 
tion between consecutive pixels. The errors in this interpolation 
can be minimized and its precision increased by selecting the 
smallest possible step. On the other hand, the computing time 
increases as the step size decreases. Therefore we settled for a 
50 m s“' step, a good compromise between computing time and 
high precision. 

The CCF width also requires particular attention because we 
require a window wide enough to cover the planet’s orbital path. 
Since the expected planet semi-amplitude is about 130 km s“^ 
(e.g., Brogi et al. 2013), we selected a window of 175 km s“^ 
on each side of the stellar radial velocity (centred on the stellar 
CCF). This allowed covering the planet’s full orbit while leaving 
on each side of the corresponding RV a continuum section large 
enough to estimate the noise level. 

3.3. Calculating the best orbital solution 

Our initial ephemerides for the orbit of 51 Peg b were taken from 
Butler et al. (2006). However, the obtained HARPS data allow 
us to derive precise radial velocities^ that can be combined with 
other measurements from the literature, so that we cover a base¬ 
line of RV measurements spanning almost 20 years with differ¬ 
ent facilities(see Table 2). 

We have thus re-derived the orbital parameters of 51 Peg b. 
These were computed using the code YORBIT (Segransan et al. 
2011). This combined dataset allowed us to derive a precise set 

^ The radial velocities derived with the HARPS DRS pipeline can be 
found in the online version. 


Table 2. Radial velocity data for 51 Peg used to derive the orbital pa¬ 
rameters. 


Instrument 

Number 

^^System 

Reference 


of points 

[km S“^] 


ELODIE@OHP 

153 

-33.2516732 

(1) 

KECK, AAT, Lick 

256 

-0.0020280 

(2) 


References. (1) Naef et al. (2004); (2) Butler et al. (2006). 


Notes. Rksystem Corresponds to the radial velocity of the system as mea¬ 
sured hy the corresponding instrument. 

Table 3. Basic orbital parameters for 51 Peg b as fitted by YORBIT 


Parameters 

Value 

f^Kystem [km S ] 

-33.152 

Period [days] 

4.231 

e 

O.O(fixed) 

a [AU] 

0.052 

htar [m S-1] 

55.2 

m 2 sin(i) [Mj„p] 

0.450 

(x> [degrees] 

O.O(fixed) 

to [BJD-2400000] 

56021.256 


of orbital parameters for the star and its planetary companion. 
The derived orbital parameters can be found in Table 3, where 
the value of RVsystem was set to the HARPS value^. During the 
fitting process, the eccentricity was fixed to zero because the ob¬ 
tained value was not statistically significant. 

3 . 4 . Recovery of the planet signal: methodology 

We extracted the planet signal from the stellar noise with the 
technique described in Martins et al. (2013). In brief, after the 
CCFs of each observation has been computed, the signal can be 
recovered in two steps: 

step 1: the CCFs are normalized with a stellar template and 
step 2: the individual CCFs resulting from step 1 are stacked 
after correcting for planetary velocity. 

Especial care in this process was taken for step 1, which con¬ 
sisted in dividing (normalizing) each of the star-i-planet CCFs by 
a carefully built stellar CCF template that represents the stellar 
signal as accurately as possible. This permitted us to remove the 
stellar signal so that only the planetary signal and noise were 
left. For this template we also needed to ensure the highest poss- 
bile S/N so that the division by this template would not introduce 
significant additional noise. To this end, we constructed two dif¬ 
ferent templates from our observations: 

Template #1 This template was constructed by stacking all the 
CCFs in our sample, after correcting each for the stellar 
radial velocity induced by the planet. This option implies 
that the template spectrum also includes a contribution of 
the planetary reflected light spectrum. Despite this fact, 
since the planet RV varies very rapidly, the planet’s minute 

^ In the fit a diiferent zero point was fitted to each instrument’s radial 
velocities. 
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Figure 1. Orbital phases of our data (black circles) and of the observa¬ 
tions used to recover the planetary signal (red stars) at maximum de¬ 
tection significance. Phase zero corresponds to the inferior conjunction 
phase, i.e., the position on the orbit where the planetary night side faces 
Earth. The green dotted line corresponds to the orbital fit from the pa¬ 
rameters in Table 3. 


signal is diluted amidst the noise and can in principle be 
neglected. 

Template #2 This template was constructed by stacking the 
CCFs of the data collected close to the expected inferior con¬ 
junction ephemerids (i.e., the position of the planet in its or¬ 
bit when its night side faces towards Earth, 0.9 < (p < Q.\, 
where (p is the orbital phase). With this template, the contam¬ 
ination of possible reflected light incoming from the planet 
is minimized. The downside of this selection is that the num¬ 
ber of available spectra for constructing the template is more 
limited than with template #1 and therefore might introduce 
non-negligible noise in the data. For comparison, since only 
20 CCFs in our sample are available for template #2, its S/N 
is ~40% of the S/N of template #1. 

To stack the CCFs after normalization by the stellar tem¬ 
plate (step 2), we discarded observations in which the planetary 
and stellar signals were spectroscopically blended or close in 
velocity. To avoid this, we only considered observations (after 
correcting for the planet’s RV computed assuming an elliptical 
orbit with the orbital parameters in Table 3) in which the ra¬ 
dial velocity difference between the planet and the star exceeds 
8 X FWHMstar (eight times the FWHM of the stellar CCF). For 
a more detailed explanation see Martins et al. (2013). Figure 1 
shows the orbital phases of our data (black circles) and of the 
observations used to recover the planetary signal (red stars) at 
maximum detection significance. 

For each individual CCF, the expected radial velocity of the 
planet was computed from the Keplerian fit presented in Table 3. 
The different processed CCFs were then re-centered by subtract¬ 
ing the planetary radial velocity from all of them (i.e., centering 
them around the velocity of the planet at any given moment) and 
co-adding them to increase their S/N. 

After the individual CCFs were stacked, a Gaussian curve 
was fitted and its amplitude was used to compute the detection 
significance. To discard fits with no physical meaning, two con¬ 
strains were set in place; 


- FWlTMpianet > 0.9 FWHMstar - This lower limit takes into 
account noise that might decrease the width of the planet’s 
CCF. For instance, if the star’s convective envelope were 
to be tidally locked to its planetary companion, the planet 
would ’’see” its host star unaffected by the star’s rotation and 
therefore the planetary CCF might be an unbroadened ver¬ 
sion of the star’s (for further details we refer to Charbonneau 
et al. 1999). This is clearly not the case for the 51 Peg sys¬ 
tem, as the stellar rotation period (~ 21 days, see Simpson 
et al. 2010) is much longer than the planetary orbital period 
(~ 4 days), and therefore this effect should be negligible. 

- FWHMpianet < 4 FWHMstar - Planetary rotation might cause 
the planetary CCF to be broadened relatively to the CCF of 
the host star while allowing for enough of a continuum for 
noise estimation. This limit was computed from the extreme 
case of a planet with twice the size of Jupiter but the same 
rotation period and observed edge-on. Nonetheless, we ex¬ 
pect the rotation period of a close-in giant planet to be much 
lower due to tidal interaction with its star. 


The significance of this detection D is then defined by 



tr noise 


( 2 ) 


where A is the amplitude of a Gaussian fit to the planet’s signal 
and (Tnoise is the continuum noise on both sides of the signal. 
We define the continuum noise as the standard deviation of the 
pixel intensity of the stacked CCFs at a separation of more than 
2 X FWHMpianet fi'om the center of the detected signal. 

As mentioned above, the planet’s real mass, and therefore 
its real orbital velocity semi-amplitude knaaeu is not known. To 
pinpoint this real semi-amplitude, we computed the detection 
significance over a range of planetary orbital semi-amplitudes 
(75-275 km s“') centered on the semi-amplitude (A:pianet ~ 133 
km s“') calculated from the minimum mass recovered from the 
radial velocity fit presented in Table 3. 


3.5. Characterizing the pianet 

When stacking the CCFs, we can expect that the maximum de¬ 
tection significance occurs when the correct velocity shift is used 
in each individual CCF that we combine. In particular, this is ex¬ 
pected to occur if we are able to input the correct semi-amplitude 
radial velocity signal of the planet. A significant detection of 
the planetary signal (and its radial velocity semi-amplitude) will 
then allow inferring the planet-star mass ratio q from 

llfpianet ^Star 

q = - = ■;-, (3) 

^Star ^Planet 

where ^pianet is the most significant amplitude as delivered from 
the analysis of the previous section, and A:star as derived in Sect. 
3.3. With q, and given the derived value for the stellar mass 
(1.04 Mq), we can compute the real mass for the planet. Together 
with the minimum mass, this value allows deriving the planetary 
orbital inclination. 


4. Results 

Following the method described in Sect. 3, we calculated the 
detection significance for evenly distributed radial velocity semi¬ 
amplitudes of the planet 75 km s“' j A:pianet i 275 km s"',, with 
a step of 0.05 km s“'. For each A:pianet we stacked the individual 
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Figure!. Detection significance as a function of A:pianet- The red line corresponds to the A:pianet value for maximum detection. The maximum 
detection occurs for similar ^pianet values for both templates. The amplitude values set to zero corresponds to ^pianet values for which no Gaussian 
fit with our restrictions could be achieved. Left panel: using template #1; right panel: using template #2. 


Planet CCF for several k2 



Radial velocity (in km/s) 

Figures. Detected signals as a function of fcpianet over the selected ve¬ 
locity range. With decreasing distance to the maximum detection, the 
signal becomes better defined and the continuum noise less dispersed. 


CCFs after correcting for the corresponding radial velocity"^. We 
then performed a simple Gaussian fit to the resulting stacked 
CCF, using the restrictions discussed in Sect. 3. The significance 
of the detection (D) was then derived using Eq. 2. This process 
was repeated for the two different template options presented in 
Sect. 3. 


The constraints mentioned in Sect. 3.4 cause the number of spectra 
used for constructing each combined CCF to vary with the assumed 
value of fepianet (it increases with kpianet)- However, the total S/N of the 
resulting CCF is always above the S/N threshold for detection of the 
observed signal (with an amplitude of 6 ± 0.4 x 10“^). 


Template #1 yields a maximum detection significance of 3.7- 
Cnoise for kpianet = 132 km s“', while with Template #2 a max¬ 
imum significance of 5.6-crnoise is computed for kpianet = 133 
km s ^ To be conservative, the error bars on the value derived 
using template #2 were computed from the 2-cr uncertainty in the 
amplitude of the Gaussian fitted to the CCF of the planet at max¬ 
imum detection significance, yielding kpianet = 

Using the same procedure to compute the error bars for the re¬ 
sults using template #1 would lead to kpianet = 

However, for this latter case, a simple visual inspection of Fig. 
2 shows that the significance curve is relatively flat for values 
between about 120 and 150 km s“* (but always above 3-(Taoise 
in this range). We thus decided to adopt all the values of kpianet 
with a significance higher than 3.0-crn(,ise, that is, kpianet = 132/;|5 
km s '. This of course lowers the significance of our detection 
to 3.0-crnoise (and not 3.7-crnoise as derived with the maximum 
value). 

Following the constraints presented in Sect. 3.4, the plane¬ 
tary CCF for maximum significance was constructed from stack¬ 
ing 25 observations of the available 90 (~27%). Although tem¬ 
plate #2 yields a higher detection significance, we decided to use 
template #1 for the remainder of the paper since it has a much 
higher S/N because many more CCFs were used to construct it. 
Therefore it will represent the stellar signal more reliably and is 
less prone to introduce additional noise into the CCF (including 
correlated noise that is difficult to quantify). 

Table 4 shows the parameters for the recovered CCF derived 
from the Gaussian fit using template #1. In this table, the ampli¬ 
tudes of the star (0.48) and planet (6 + 0.4 x 10“^) CCFs coiTe- 
spond to the depth of the CCF when the the continuum has been 
normalized to 1. The planet-to-star flux ratio will then be given 
by the ratio of these two quantities. The significance of the de¬ 
tection corresponds to the maximum significance for 75 km s * j 
^Planet i 275 km s“* and was computed by dividing the amplitude 
of the planetary CCF by the noise on the wings of the CCF. The 
FWHM is the full width at half maximum of the fitted Gaussian, 
7.43 km s“' for the star, 22.6 + 3.6 km s“' for the planet. 

Figure 3 shows for several assumed values of kpianet over the 
adopted range a section of the stacked CCFs centered on the ra¬ 
dial velocity where the planetary signal would be expected from 
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Planet CCF for several k2 



-32-24-1G-8 0 8 16 24 32 -32-24-16-8 0 8 16 24 32 -32-24-16-8 0 8 16 24 32 

Radial velocity (in km/s) 


Figure 4. Detected signals as a function of fepianet > but for random sec¬ 
tions of the normalized CCF where the planetary signal is not expected 
to be found. 


a Keplerian fit to the observations using the parameters in Table 
3 and the selected A:pianet- The plot illustrates how the recov¬ 
ered CCF changes when we co-add spectra for different assumed 
semi-amplitude values. When a value for the semi-amplitude 
close to 132 km s“' is assumed, the recovered signal is well 
defined and its wings are less noisy; this is closer to the distribu¬ 
tion of the expected Gaussian shape of a CCF. However, as the 
assumed value of the semi-amplitude departs from 132 km s ', 
the CCFs of the individual observations will be co-added to each 
other, albeit imperfectly. Therefore, in these cases a signal is also 
expected to be detected, but spread out across the continuum and 
of lower signihcance because the wings additionally contribute 
to noise. When close enough to the best-ht value of the semi¬ 
amplitude, a signal above the noise level can still be seen. This 
signal might seem to be of similar amplitude as the one detected 
for the correct value of orbital semi-amplitude, but it will have a 
lower signihcance due to the increased noise in the wings. This 
can be seen in the panels for A:pianet = 115 and 120 km s“', which 
show signals that might seem to be of similar amplitude to the 
one in the central panel, but because of the increased noise in the 
wings, these signals are of lower signihcance according to Eq. 2 
(see Fig. 2). 

To test whether the observed signal might be a spurious com¬ 
bination of random noise, we repeated the process, but instead 
of using a Keplerian function to compute the expected radial ve¬ 
locity signal of the planet, we attributed a random radial velocity 
within the range [-50, 80] km s * to each CCF. The result of that 
analysis can be seen in Fig. 4. They show that no signihcant sig¬ 
nal is detected, as expected. This test gives us conhdence that the 
detected CCF of the planet is not a mere artifact. The dips that 
appear in the panels of Fig. 4 appear at random positions, and 
their FWHM is always lower than the FWHM of the star. For 
instance, the FWHM of the center dip in the 115 km s ' panel of 
Fig. 4 has a FWHM of only 3.8 km s“'. This strongly suggests 
that these peaks are purely caused by noise, especially as they 
appear randomly positioned. The constraints that we have speci- 



FigureS. Planet signal for maximum detection and fitted Gaussian 
curve with the Levenberg-Marquardt method. 

fied in Sect. 3 ensure these dips as discarded as nonphysical and 
are not confused with the planetary CCF. 

To verify whether our analysis data-processing was sound, 
we simulated 100 sets of noiseless CCF functions simulating 
idealized observations of the star -H planet signal for random val¬ 
ues of the planetary semi-amplitude A:pianet in the range [100, 
180] km s“'. These sets of simulated observations consisted 
of 90 star-tplanet CCFs (the number of observations we have), 
where the planet radial velocity is computed from one of the ran¬ 
domly selected A:pianet through Eq. 3. Each simulated CCF was 
built by positioning the stellar Gaussian at the observed velocity 
of the star at each given moment (obtained from the real CCFs) 
and co-adding the planet Gaussian with an amplitude of 5 x 10“^ 
(a value similar to the expected for our test subject) . The goal 
was to verify whether the injected signal was always fully recov¬ 
ered by the reduction process. The results show that the injected 
signal was successfully recovered in all simulations. The ob¬ 
tained values of A:pianet were always the values that were injected 
into the expected figure, with a standard deviation of only 0.11%. 
Similar results were obtained for the amplitude, which was al¬ 
ways recovered with a standard deviation lower than 0.01%. This 
test shows that that the data-analysis pipeline works correctly 
and can be used safely to retrieve the planetary signal. 

The planetary CCF with the highest recovered significance is 
plotted in Fig. 5 together with a Gaussian fit using a Levenberg- 
Marquardt algorithm. After removing the planetary CCF and fit¬ 
ting it with a Gaussian, we used the residuals of the fit to de¬ 
fine the error bars. To do so, we injected into the residuals a 
Gaussian curve with the same parameters as the detected plane¬ 
tary CCF after subtracting the detected Gaussian, but at different 
radial velocities. This procedure was repeated 10000 times, and 
in each we recovered the signal injected by fitting a Gaussian 
profile. The standard deviation of the recovered FWHM and am¬ 
plitude values were then considered to be the 1-cr errors as listed 
in Table 4. 


5. Discussion 

Our results suggest that we were able to successfully detect the 
light spectrum of 51 Peg reflected on its hot Jupiter companion. 
The results also indicate that the best-fit semi-amplitude of the 
orbital motion of the planet is A:pianet = 132]^[^ km s“'.. 
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Table 4. Comparison of the stellar and planetary CCF parameters. 


Parameter 

Star 

Planet 

Amplitude 

0.48 

6.0+0.4 X 10^5 

Significance [cTnoisJ 

- 

3.7+0.2 

FWHM [km s^'] 

7.43 

22.6+3.6 


Notes. Comparison of the stellar and planetary CCF parameters. For the 
star signal, we present the median value of the amplitude and FWHM 
of its CCF for all observations. For the planetary CCF, we present the 
values of the amplitude and FWHM of its Gaussian fit and its detection 
significance. In both cases, the level of the CCF continuum flux has 
been set to one. 


The planetary mass cannot be lower than its real mass (corre¬ 
sponding to an orbital inclination of 90 °), which places a higher 
limit on the planetary orbital semi-amplitude of A:pianet ~ 133 
km s"'. Combining the detected value of the semi-amplitude 
with this assumption, Eq. 3 yields a real mass for 51 Peg b of 
0.46;^° qj Mjup for a stellar mass of 1.04 Mq (Santos et al. 2013). 
By comparing this with the derived minimum mass of m 2 sin i = 
0.45 Mjup (Table 3), we can infer an orbital inclination of 80)^}® 
degrees. This result is compatible with the results obtained inde¬ 
pendently by Brogi et al. (2013) - 79.6° < I < 82.2° for a 1 cr 
confidence level. 

These results suggest that we were able to successfully re¬ 
cover the reflected planetary signature, but the noise present in 
the data hinders constraining parameters like the planet’s geo¬ 
metric albedo. It is interesting, however, to verify which possible 
values of the albedo we would recover if we assumed plausible 
values for the radii of the planet according to Eq. 1 and assumed 
that the planet-star flux ratio is given by the ratio of the ampli¬ 
tudes of the the detected planet and stellar CCEs. Note that the 
scatter of radii observed for hot-Jupiter planets of a given mass 
does not allow us to simply use the derived planetary mass to 
estimate a reliable value for the planetary radius (e.g., Baraffe 
et al. 2010; Eortney & Nettelmann 2010). 

The results show that the obtained Ag is above unity for a 
radius of 1.2 Rjup . Eurthermore, if the planet’s reflectivity is 
assumed to be Lambertian (i.e., the planetary surface reflects 
isotropically in all directions), it has been shown that the plan¬ 
etary albedo will have at most a value of two thirds (Sobolev 
1975). Higher values imply that the planetary surface or atmo¬ 
sphere would be strongly backscattering. It has been shown that 
cloudy planets in the solar system tend to be close to Lambertian 
(e.g., Venus - see Traub & Oppenheimer 2010), while rocky 
planets tend to be more uniformly bright than a Lambert’s 
sphere (e.g., Earth and Moon - see Traub & Oppenheimer 2010). 
Moreover, literature values suggest that hot Jupiters usually have 
albedos lower than ~0.3 (Cowan & Agol 2011). However, sev¬ 
eral works (e.g., the case of HD 189733b, Evans et al. 2013; Pont 
et al. 2013) repotted albedo values as high as 0.4-0.5. Such high 
values may be explained by the scattering of condensates in the 
high atmosphere (Sudarsky et al. 2003). 

Eor our measurements, the observations used to recover 
the planetary signal were obtained at an average orbital phase 
g(a) = 0.87. Thus, for a perfectly reflecting Lambert sphere, 
that is, Ag = 2/3, Eq. 1 yields a radius for the planet of 1.6+0.2 
Rjup (the error bar denotes the error found in the detected CCE 
amplitude). Eor an albedo of 0.5, the planetary radius will in¬ 
crease to 1.9+0.3 Rjup. Such high radii have been observed in 


other hot Jupiters of similar mass and orbital period as 51 Pegb 
(e.g., Kepler-12b, Eortney et al. 2011), making our assumptions 
plausible. 

The estimates discussed here are based on a 3-sigma detec¬ 
tion, however. As we can see in Eigure 5, the detection we ob¬ 
tained shows a strongly correlated noise signature, which is dif¬ 
ficult to take into account when fitting the observed signal with 
a Gaussian function. We can thus not fully trust the derived pa¬ 
rameters for the detected planetary CCE. This is crucial for the 
EWHM of the recovered CCE, which might be used to estimate 
the planetary rotation velocity (e.g., Santos et al. 2002). Such a 
high CCE width, when compared with the stellar CCE, would 
imply an extremely rapid rotation for the planet (~18 km s '), 
much higher than expected for a planet assumed to be tidally 
locked to its host star (~ 2 km s“'). Nonetheless, because of 
our low confidence in the recovered CCE parameters, any con¬ 
clusion regarding the planetary rotation velocity would be too 
speculative. 

An alternative explanation for the broadening might be an 
imperfect orbital solution. But our target has been followed for 
many years, and the star’s orbital parameters are very well con¬ 
strained. Eurthermore, the errors in the stellar RV are on the or¬ 
der of the ms*, which transposed into the planet RV domain 
would yield at most errors in RV of 2-3 km s *. This is insuffi¬ 
cient to explain such broadening. 

Even though there is no stellar (visual) companion to 51 Peg 
reported in the literature (e.g., see Luhman & Jayawardhana 
2002), we explored the possibility that the observed signal might 
be caused by the spectrum of a contaminating background star. 
We used the Besangon model (Robin et al. 2003) to determine 
the probability of a stellar contaminant close to 51 Peg, whose 
spectrum would cause the magnitude of the observed signal. 
Given the area of the CCE for the (here suggested) planetary sig¬ 
nal (10 "*) and the magnitude of 51 Peg, the contaminant should 
be of magnitude 16 or brighter (brighter if the contaminant does 
not fall entirely within the HARPS fibre, which has a diameter of 
1 arcsec). According to the Besangon model, there are only 760 
EGK and M stars per square degree at the position of 51 Peg b 
that possess a magnitude brighter than this value^ . If we conser¬ 
vatively assume that any contaminant at a distance smaller than 
2.6 arcsec were able to produce the observed signal (for details 
see Cunha et al. 2013), we can expect a probability of a stellar 
contaminant of at most only 0.1%. We can thus reasonably con¬ 
fidently discard the possibility that the signal is produced by a 
background stellar contaminant. 

Note that the sky brightness has a magnitude per arcsec^, 
in V band, of 21.8 on a new moon night, and of 20.0 on a full 
moon night (e.g., see Cunha et al. 2013), which excludes this as 
a possible source of contamination. 

A small offset can also be seen on the center of the CCE 
recovered with the highest detection significance, which can be 
expected to be centered around zero for the real kpianet- Several 
explanations can be found for this offset: i) it might be phys¬ 
ically induced, for instance by high-speed winds from the day 
to the night side in the planetary atmosphere ( see Miller-Ricci 
Kempton & Rauscher 2012); ii) it might be an observational bias 
that results from stacking observations of the planet all on the 
same side of the superior conjunction while not at full phase (see 
Eig. 1); iii) and it might be a higher contribution of noise in one 


^ We only chose FGK and M stars since other spectral types would 
not produce a CCF detectable by the HARPS reduction, which used a 
G2 template for the cross-correlation. 
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RV (km/s) 

Figure 6. Simulated signal injected at +60kms“* in relation to the po¬ 
sition of the center of the planetary CCF. 

of the sides of the CCF. However, we consider that this offset is 
not significant because of the high noise. 

To test how much the observed signal can be affected by this 
noise, we performed the following test. For each observed CCF 
we injected an artificial signal (a pure Gaussian) with an ampli¬ 
tude of 5 X 10“^ (similar to the retrieved planetary signal), but 
at a different position in velocity (at H-60kms * when compared 
to the planet). We then ran our data-reduction process and tried 
to retrieve this signal. As can be seen in Fig. 6, we recovered the 
injected CCF with a similar amplitude of 4.6 x 10“^. However, 
the FWHM of the fitted Gaussian is ~ 27.6 km s ', significantly 
larger than the injected FWHM value (7.43 km s'* as observed 
in the stellar CCF). This test was repeated for different signal 
intensities, and we found that when the signal strength is much 
higher than the noise, it is recovered with parameters close to the 
injected values. Nonetheless, when the signal strength is close to 
the noise level, the recovered signal appears broadened. This test 
shows that the parameters of the recovered CCF are strongly af¬ 
fected by the noise. 

In this particular simulation we show that the obtained 
FWHM is much broader than the injected one. Interestingly, the 
value for the FWHM we recovered in our data for the plane¬ 
tary signal (Table 4) is also much higher than the value observed 
for the stellar CCF. This shows that the FWHM we obtain can¬ 
not be considered reliable to derive the rotational velocity of the 
planet or the velocity of its atmospheric winds, for instance (e.g., 
Snellen et al. 2010). 

6. Conclusions 

We presented the first application of the technique described by 
Martins et al. (2013). We were able to find evidence for the re¬ 
flected light from the hot-Jupiter 51 Pegb with this technique. 
This result is encouraging and constitutes a very valuable proof 
of concept. Our method can be used to recover an exoplanet’s 
spectroscopic reflected signature from among the stellar noise, 
despite the extremely low planet-to-star flux ratio. 

Although only a 3 cr„oise detection was possible, the ob¬ 
tained data allow us to place constraints on the planetary mass 
(0.46j;QQj Mjup) and orbital inclination (80^}g degrees), assum¬ 
ing the detection is real. As presented above, reasonable values 
for the albedo (0.5) are found if we assume that 51 Peg b is 


an inflated hot Jupiter with a radius of about 1.9 Rjup. Strongly 
inflated planets with similar mass and orbital period have been 
reported in the literature (with radii measured with transits). 

Of particular interest is the fact that this detection was possi¬ 
ble with data collected with a currently existing observing facil¬ 
ity on a 4-meter-class telescope (the HARPS spectrometer at the 
3.6 m ESO telescope). Additional data will also allow addition¬ 
ally constraining of some of the inferred parameters by increas¬ 
ing the detection significance and allowing for a better fit to the 
detected planetary CCF. 

These encouraging results clearly show a bright future for 
this type of studies when next-generation instruments (e.g., 
ESPRESSO at the VET) and telescopes (e.g., ESO’s E-ELT) be¬ 
come available to the community. The sheer increase in preci¬ 
sion and collecting power will allow for the detection of reflected 
light from smaller planets, planets on orbits with longer periods, 
or an increase in detail for larger planets like 51 Peg b. 
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